A singular Hamilton-Jacobi equation modeling the tail problem 
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Abstract 



^ Wc study the long-time/long-range behavior of reaction diffusion equations with negative square 

root-type reaction terms. In particular we investigate the exponential behavior of the solutions after 
space and time are scaled in a hyperbolic way by a small parameter. This leads to a new type 

rG ' of quasi-variational inequality for a Hamilton-Jacobi equation. The novelty is that the obstacle, 

which defines the open set where the solutions of the reaction diffusion equation do not vanish in 
the limit, depends on the solution itself. Counter-examples show a nontrivial lack of uniqueness for 
the variational inequality depending on the conditions imposed on the free boundary of this open 
set. Both Dirichlet and state constraints boundary conditions play a role. When the competition 
term does not change sign, we can identify the limit while, in general, we only obtain lower and 

£vi ■ upper bounds. 

Although models of this type are rather old and extinction phenomena are as important as 
blow-up, our motivation comes from the so-called "tail problem" in population biology. One way 
to avoid meaningless exponential tails is to impose a singular mortality rate below a given survival 

t^- ■ threshold. Our study shows that the precise form of this singular mortality term is asymptotically 

irrelevant and that, in the survival zone, the population profile is impacted by the survival threshold 
(except in the very particular case when the competition term is nonpositive). 
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1 Introduction 

We study the asymptotic behavior, as e — > 0, of the solutions to reaction-diffusion equations (with 
singular reaction term) of the form 



- eAn £ = \n e R - ±{p £ n e ) l/2 in R d x (0, +oo) 
e «2/e on R d x {0}, 



n E,t /i\ 
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where u® : R rf — > R is a given function, R : M. d — > R represents a linear logistic growth/death rate and 
the survival threshold parameter /3, which models a singular death term, is given, for some u m < 0, 
by 

p e = e Um ' e . (2) 

The positive parameter e is introduced by a hyperbolic scaling (x, t) t— > (x/e, t/e) with the aim to 
describe the long time and long range behavior of the unsealed problem (corresponding to e = 1). 
The limiting behavior of scaled reaction-diffusion equations with KPP-type reaction has been studied 
extensively in, among other places, the theory of front propagation (U [T9| I12j ) using the so called 
WKB- (exponential) change of the unknown. 

The novelty of the problem we are considering here is the presence of the negative square root 
term. To the best of our knowledge, the first study of such nonlinearity goes back to |11| where it 
is proved that local extinction occurs, i.e., the solution can vanish in a domain and stay positive in 
another region. For this reason (3 is thought to represent a survival threshold. That a solution of a 
parabolic problem can vanish locally is a surprising effect and as singular as the blow-up phenomena 
for supercritical reactions terms ([IZ])- In population biology such behavior prevents the so-called 
"tail problem" where very small (and thus meaningless) populations can generate artifacts (|14j). 
Although the mathematical analysis of the limit of ([1]) turns out to be a full subject in itself, our 
primary motivation comes from qualitative questions in population dynamics. 

Indeed ([I]) is the simplest model for studying the effect of "cutting the tail" but many other problems 
are relevant in ecology. Along the same lines, in the context of front propagation, one may consider 
the modified Fisher-KPP equation 

n E ,t - eAn £ = -n e (l - n £ ) (/3 £ n £ ) 1/2 in M. d x (0, +oo), 

and ask the question whether the square root term changes fundamentally the study in |12] and [H] of 
the propagation of the invading/combustion fronts. In the context of speciation, an elementary model 
in adaptive evolution is the non-local reaction-diffusion equation 



n £ ,t 



11 [ 

eAn £ = -n £ R(x,I £ ) ((3 £ n £ ) 1/2 in R d x (0,+co) with I £ (t) = / ip(x)n E (x,t)dx, 



where n £ is the population density of individuals with phenotypical trait x, R represents the net 
growth rate, ip is the consumption rate of individuals and I(t) is the total consumption of the resource 
at time t. The survival threshold was introduced in p3]. Finally e may represent large time and small 
mutations as studied in [H EJ [16] . It is known that under some assumptions the density concentrates as 
an evolving Dirac mass for the fittest trait. In biological terms this means that one or several dominant 
traits survive while others become extinct. Phenomena such as the discontinuous jumps of the fittest 
trait, non smooth branching and fast dynamics compared to stochastic simulations, motivated [14] to 
improve the model by including a survival threshold. Numerical results confirm that this modification 
gives dynamics comparable to stochastic models. It is interesting to investigate rigorously whether 
the dynamics of the Dirac concentration points are really changed by the survival threshold and to 
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explain why its specific form (n e versus nj with < 7 < 1) seems irrelevant. 

A way to approach these questions for ([T]) is through the asymptotic analysis of n £ . Since, as in 
the classical case, i.e., the Fisher-KPP equation without the square root term (see JT2]), n £ decays 
exponentially, the limit is better described using the Hopf-Cole transformation 

u £ = elnn £ , (3) 



which, for u® = elnn®, leads to the "viscous" Hamilton- Jacob i initial value problem 



u e ,t-eAue-\Du e \* = R-expl "^ £ j in R d x (0, +00), 
u° in M d x {0}. 

Throughout the paper we assume that there exist C > and u° S C 0,1 (1R ) such that 

||i?|| c o,i < C and ||u°||co,i < C, (5) 

u° £ £ C(M d ), n°<C and u° — ► u° in C(M. d ). (6) 

In the limit e — > 0, it is easy to see, at least formally, that any local uniform limit of the family 
(u e ) e> o will satisfy, in the sense of the Crandall-Lions viscosity solutions (|10|). the Hamilton-Jacobi 
free boundary problem 

'u t = \Du\ 2 + R in nc R d x (0,oo), 
u = -00 in ff n (R d x (0, 00)), 
u > u m in 0, 
u = n° in Hn(IR d x {0}), 

with the space-time open set fi defined by 

Q =lnt {(x,t) £ R d x (0,oo) : lim u e (x,t) > -00}. 

Notice that ([7]) is an obstacle problem with an obstacle depending on the solution itself. As a 
matter of fact the open set O plays an important role and, hence, the problem may be better stated 
in terms of the pair (u,Q). The difficulty is that ([7]) has several viscosity solutions (see Appendix lAl 
for examples) depending on the sense the boundary conditions are achieved and the sign of R. 

Next we discuss the two boundary conditions arising in d7J). The first is the Dirichlet boundary 
condition in the third relation in (J7J). Its precise form is 

lim u(x,t) = u m . (8) 

(x,t)-t(xo,to)£dSl 

The second is the state constraint boundary condition (see [IS]), which is natural in view of the second 
equality in ([7]). It states that 

u is a supersolution in f2 and a subsolution in £1 . (9) 

The basic questions we are considering in this paper are: 

• What boundary condition should be satisfied by the limits of the family (u e ) e> q on <9f2? Dirichlet or 
state constraint? The latter appears to play a fundamental role. To the best of our knowledge, there 
are no results available for state constraint problems with time varying and non smooth domains. 
Most of the technicalities in the paper stem from this difficulty. 

• Does the limit e — > select a particular solution to (|7|), i.e., is there a natural selection? Is the limit 
of the family (u e ) e >o the maximal subsolution or minimal solution to dTJ)? 

• Do the limits of the family (u e ) e> Q depend on the specific form of the survival threshold, i.e., can 
we replace (/Seng) 1 ' 2 by (/3 £ n £ ) 7 with 7 € (0, 1) without affecting the outcome? 



An important ingredient of our analysis is the asymptotics, as e — > 0, of the solution u\ of 

in M rf x(0,+oo), 
x{0}, 
which is obtained, after the Hopf-Cole transformation 



u\ >t = £& u l + \ Du l\ 2 + R 



u\ = u° £ in ™ d 



u £ = e In n £ , 



from the simplified reaction diffusion equation 



nit 



>i: 



eAnl 



-^R in 



x (0, +00), 



(10) 



(11) 



(12) 



exp(e -1 ^) in R d x {0}. 



In view of flSJ) and (|SJ), it follows from [12] that, as e — > 0, the sequence («g) e converge locally 
uniformly to u 1 S C(M d x (0, 00)), which is the unique viscosity solution of the eikonal -type equation 



} = \Du 1 \ 2 + R in R d x(0,+oo), 



it 



u 1 = u° in 



(13) 



x{0}. 



The maximum principle yields n £ < n\, which in turn implies that u £ < u\ and, in the limit (this is 
made precise later), u < u . It also follows from (j4]), at least formally, that, as e — > 0, 



-00 



m 



x (0,oo))\^ 1 , 



where 



n l 



{(x,t) I v}(x,t) > u m }. 



(14) 



It turns out that the case of nonpositive rate R is particularly illuminating and the above ques- 
tions can be answered completely and positively using u 1 (see Section [2]). The problem is, however, 
considerably more complicated when R takes positive values. In this case we introduce an iterative 
procedure that builds sequences of sub and supersolutions (Section [3]). This construction gives the 
complete limit of u e when R is constant (Section H]). The limit is not the maximal subsolution of (|7|) 
and the Dirichlet condition is not enough to select it. In Section \5\ we consider strictly positive spa- 
tially dependent R and provide a complete answer in terms of the iterative procedure. The relative 
roles of the Dirichlet and state constraint boundary conditions appear clearly in this case. In Section 
[6] we summarize our results. In the three part Appendix we present some examples of nonuniqueness 
as well as the proofs of few technical facts used earlier. 

We conclude the introduction with the definition and the notation of the half-relaxed limits that 
we will be using throughout the paper. To this end, if (w e ) e> o is a family of bounded functions, the 
upper and lower limits, which are denoted by w and w respectively, are given by 



w(x) = lim sup w; £ (y) 
e— >0,y— >x 



and w(x) 



liminf w e (y). 



(15) 
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2 Nonpositive growth rate 

Here we assume 

R < in R d , (16) 

and show that the behavior of the family (u e ) e , in the limit e — > 0, can be described completely in 
terms of solution u 1 of fjl3j) . which carries all the necessary information. More precisely, we can state 
the 

Theorem 2.1. Assume {3Jj ; (0|) and dip]). ^4s e — > 0, t/ie family (u E ) e> o converges, locally uniformly 
in n 1 and in (R d x (0,oo)) \W, to 

u(xt)= f ul &t) for (x,t)en\ 

\ -oo /or (x,t) G (M d x (0,oo))\O 1 , 

wi/i u 1 and il 1 defined by Iil3\) and fli^| ) respectively. In particular, u(x,t) — > u m as (x,t) — > ^fi 1 . 

Before we begin with the proof, we present and discuss below several remarks and observations 
which are important to explain the meaning of the results. 

Firstly, by "uniform convergence" to — oo, we mean limsup e ^ 0y _> xs _ >t u £ (y,s) = — oo. Secondly, 
the u associated with the open set Q 1 is the maximal solution to ([7]). Indeed any other solution u, 
with the corresponding open set Q, satisfies u < u l and thus £1 C SI 1 and u < u. It also satisfies the 
Dirichlet and state constraint boundary conditions. To verify the latter we notice, using the standard 
optimal control formula ( [15} [T3| [T] ) , that 

u 1 {x,t)= sup { I ( -!?W- + R( x (a)) J ds + u (x(0)) : x €C 1 ([0,t];R d ) 

If x(-) is an optimal trajectory, the dynamic programming principle implies that, for any < r < t, 

U \x,t)= [ l-^^- + R(x{s)))d S + U 1 (x(T),T). 



Since R is nonpositive, u 1 is decreasing along the optimal trajectory. It follows that, if n 1 (x, t) > u r 
then, for all < r < t, n 1 (x(r),r) > u m . 
Hence, for all (x,t) G O , 

u(x,t)= sup {f(-—^- + R(x(s)))ds + u (x(0))-.x£C 1 ([0,t]-,W 
(itsj.sje!! 1 Uo V 4 / 

x(t) — X 

and, therefore, u verifies the state constraint condition. 
Finally, the limit u does not depend on the details of the singular death term. In particular it is 

1 /9 

the same if we replace in ^Q) n e exp(e _1 n m ) by njexp(e _1 7n m ) with < 7 < 1. Hence, the value 
7 = 1/2 is irrelevant. 

We continue with the 



Proof of Theorem \2.1\ As already discussed in the introduction, we know that u £ < u\ but we cannot 
obtain directly the other inequality in the limit e — > 0. It is therefore necessary to introduce a pair of 
auxiliary functions v A and v £ ' which converge, as e — > 0, in C(M. d x (0,oo)) to max(u 1 , —A). Using 
this information for appropriate values of the parameter A, we then prove that, as e — > 0, u £ — > u l 
locally uniformly in the open set 

A = {(x,t) : u 1 ^,*) >u m }, (18) 

and u £ — > — oo locally uniformly in the open set 

B = {{x,t) : u\x,t) <u m }. (19) 

To this end, for any A such that 

< A < -Um, (20) 

we consider the functions v A and v £ ' given by 

n £ + exp( ) = exp( — ) and n £ + exp( ) = exp( ). (21) 

We have: 

Proposition 2.2. Assume $5j), ^j), [To}) and [20]) . As e — > 0, the families (v £ ' ) £ >o and (v^) £> o 
converge in C(R x [0,oo)) to the unique solution v ' = max(n 1 ,— A) of 



mm [v 



AA , A „, AA 



+ A,Vt>-\Dv A < 1 \ 2 -R)=0 in M d x(0,oo) 



(22) 
) AA = max(u°, -A) on R d x {0}. 

We postpone the proof to the end of this section and next we prove the convergence of the family 
(u £ ) £ in the sets A and B. We begin with the former. 

Fix (xo,to) E A. By the definition of A we have n 1 (xo,to) > u m and, hence, we can choose A such 
that u 1 (xo,io) > ~A > u m . Proposition 12.21 yields that, as e — > and uniformly in any neighborhood 
of (x ,t ), 



, A -> v A ' 1 



max(— A,u ) = u . 



Using the latter, the choice of A and the fact that 

u £ = v £ +eh\{\ -exp(e _1 (-,4- v A ))), 

we deduce that, as e — > 0, u £ — > u 1 uniformly in any neighborhood of (xq, to). 

Next we consider the limiting behavior in the set B. To this end, observe that, using ([3]) and (llip . 
we find u £ < ul and, thus, passing to the limit in the viscosity sense, u < u 1 and 

u < u m in B. 

Assume that, for some (xo,to) G B, u(xo,to) > — oo. Since u is upper semicontinuous (see [2]), there 
exists a family (4> a )a>o of smooth functions such that u — (j) a attains a strict local maximum at some 
(x a , t a ) and, as a — > 0, 

(x a , t a ) ->■ (x , to), u{x ai t a ) > u(x ,t ) and u(x a , t a ) -)• u(x , t ). 



It follows that there exists points (x Qj£ , t a , s ) such that u £ — <j) a attains a local maximum at (x a ^ e ,t a ^ e ), 
(x a ,e,t a ,e) -> (»a, t a ) as e -> 0, and, in view of (|4|), at (x a>e ,t a>e ), 

<f) att -sA(f) a - \Dcj) a \ 2 - R < -exp((2e) _1 (u m - u £ )). 

Letting e — >• we find that at (x Q , i a ) 

a ,t - \D(j) a \ 2 - R < limsup[-exp((2e) _1 (n m - u £ (x at£ ,t a>£ )))]- 

The definition of u yields 

limsupu £ (x ai£ ,t Q)£ ) <u(x a ,t a ) 

and, since, for a sufficiently small, u(x,t) < u m , we have 

u(x a ,t a ) < u m and limsup[- exp((2e) _1 (u m - u £ (x a>£ , t a>£ )))] = -oo 

and, finally, at (x a ,t a ), 

(j> a ,t - \D(j) a \ 2 - R < -oo, 

which is not possible because <j) a is a smooth function. 

The claim about the uniform convergence on compact subsets is an immediate consequence of the 
upper semicontinuity of u and the previous argument. □ 

We conclude the section with the proof of Proposition 12.21 Since it is long, before entering in the 
details, we briefly describe the main steps. We begin by establishing independent of e bounds on 
the family (v£) e . Then we show that the half-relaxed limits v a and v a are respectively sub and 
supersolutions of (|22p . We conclude by identifying the limit. 

Proof of Proposition [27B . By the definition of v £ , we have v £ > —A and, thus, the family (v £ ) £ is 
bounded from below. 
To prove an upper bound we first notice that, on M. d x {0}, 

A n — A— 1»2 . A+u® 

vf = u o + em ( 1 + e — ^) and v f = -A + e\n{l + e-^), (23) 

hence, 

i^<max(it° + eln(2),-A + eln(2)) on R d x {0}, 

and, finally, in view of ([6]), 

vf < C A on R d x {0}, 



for Ca > such that max(— A, u £ ) < Ca- 
Moreover, since R < 0, we have 



n e „ (Psn £ ) 1/2 

— A — A 

n £ + exp(— ) n £ + exp(— ) 



v? t - eAvf - \Dv?\ 2 = ,-a. R ,-a^ ^ ° in M x (°» °°)' ( 24 ) 



It follows from the maximum principle that 

v £ < C A + eln(2) in R d x(0,oo) 



Next we show that v is a supersolution of (|22p . Since u m < —A and 

(Pen e ) 1/2 . (Z^) 1 / 2 1 ,n m + A, 

^ TT7? = o ex P(" 



^ + exp(^)- 2ri£exp( ^i ) 1 /2 2 2, 

as e — )• and uniformly on M. d x (0, oo), we have 

{Pen £ ) 1 ' 2 



ra £ + exp( — ) 



0. (25) 



From dUD, dMD and 



^ ^T < I, 

n £ + exp(- 

we then deduce that, in R d x (0, oo), 



v A t - eAv A - \Dv A \ 2 >R- O(e), (26) 

while by the definition of v A we also have 

v A + A>0. (27) 

Combining (J26|) and ((271) and using the basic stability properties of the viscosity solutions (see [2]) 
we find that the lower semicontinuous function v is a viscosity supersolution of (1221) . 

To prove that v is a subsolution to (|22p , following classical arguments from the theory of viscosity 
solutions (see [2]), we fix a smooth (/> and assume that v — </> has a strict local maximum at (xo,to). 
It follows that there exists a family, which for notational simplicity we denote again by e, of points 
(x £ ,£ £ ) £ >o in M. d x (0, oo) such that v A — <j) has a local maximum at (x £ ,£ £ ), and, as e — >• 0, (x £ ,t £ ) — >• 
(x ,to) and v A (x £ ,t £ ) -^v A (x ,t ). 

We also know, still using (j2l|) and (j2"5|) . that v A solves 



.4 



v 



A t - eAttf - |^| 2 = (1 - exp(^-^)) J R - 0( e ) 



It then follows that, at (x £ ,t £ ), 

<f>t - eAcj) - \Dcj)\ 2 - (1 - exp(e- 1 (-^ - v A ))R < O(e). (28) 

Recall that lim £ _>ov A (x £ ,t £ ) = v A (xo,to) > —A. Hence, ifv A (xo,to) > —A, then 

]imexp(e- 1 (-A-v A (x £ ,t £ ))) =0. 

From this and (J28|) we deduce that, ifv(xo,to) > —A, then, at (xo,£o)> 

4> t - \D<f)\ 2 -R<0. 

Next we show that v and v satisfy the appropriate initial conditions. Indeed, in view of ([6]) and 
23|) . we know that, as e — > 0, 

^ ->• max(-^,u°) on M d x {0}. 



It also follows from a classical argument in theory of viscosity solutions ([Sid]) that, on M. d x {0}, 

v A - max(-A, u°) < and v A - max(-A, u°) > 0. 

and, hence, v and v satisfy respectively the discontinuous viscosity subsolution and supersolution 
initial condition corresponding to (|22p . 

We already know from the definition of v and v that v < v A , while from the comparison property 
for (l22j) in the class of semicontinuous viscosity solutions (see [TJ EJ [9] ) we conclude from the steps 
above that v < v in M. d x (0, oo). Hence v = v = v ■ is the unique continuous viscosity solution 
of (|22p and, consequently, the families v A and v e ' converge, as e — > and locally uniformly, to v ' 1 . 

Combining ([3]) and (|2ip we find 



V ^ L =u L £ +eln(l + exp^ ) and vf' v = -A + eln(l + exp^r^). 

Moreover, as we already explained it in the introduction (see (|12p - (|13p ). we know that, as e — > 0, 
u\ — > u 1 locally uniformly. Hence, always for A < —u m , we obtain that, as e — > 0, 



v 



£ 



max(u , — A) locally uniformly in K x [0, oo). 



It also follows that the family (v a ) £> q converges, as e — > 0, locally uniformly to v > = max(u 1 , —A). 

D 

3 General rate 

When R changes sign, the situation is much more complicated and (|17p does not hold in general. 
In this case we are able to provide only inequalities for the half-relaxed limits of the family (u £ ) £> o- 
These estimates are used later to characterize the limit when R is positive. 

Fix uq, 5 > and recall that u l is the solution of (|13p with u 1 = uq on M. d x {0}. We introduce next 
the family (uf [uo],Cf [uo],^f[uo])i e z+ which is defined iteratively. To this end, for i = 1, let 

u*[u ] = u 1 , C([u ] =^ d x [0,oo) and fif [u ] = {(x,t) e R d X [0,oo) : u{ [u ](x,t) > u m - 5}, (29) 

and, given uf [uo],Cf [uq] and £lf{uo], uf +1 [uo] : R d x [0, oo) — > R U {— oo} is defined by 

|x(s)| 2 



n f+iM(M) = sup I j 



+ R(x(s)) 



ds + Uq (x(0)) : 



(30) 



x G C 1 ([0,t];M d ), (x(s),s) G f^[u ] for all s G [0,t], x(t) = x} 
with 



CffiN] = {(x,<) G nf[«o] : uf +1 [« ](x,t) > -oo} (31) 

and 

^+iM = {(x,t) G nfluo] : uf +1 [u ](x,t) >u m -5}c Cf +1 [u }. (32) 

It follows that, in general, Cf +1 {uo] C il^[uo]. The inclusion may be, however, strict, i.e., they may 
exist points (x,t) G ^f[«o] which cannot be connected to M x {0} by a C 1 trajectory staying, for all 
s G [0,t] in fif[u ]. (See Figure H) 



n?[«o] 




Figure 1: An example of the space-time set f^[iio]. The point (x,t) G ^f[^°] cannot be connected to 
R rf x {0} by a C 1 trajectory (x(s), s) s ^[o,t] staying within fi|[tt ]. More generally, for the points in the 
grey area, called V, there is no admissible trajectory. We have indeed Cf +1 [i*o] = &f [^o] \ U- 



Moreover jSJ) , (j29|) and classical considerations from the optimal control theory ( \15\ [T3"l [TJ [8] ) yield 
that, for all i G Z + , the sets Cf [u$] and flf are open and uf [uq] G C(C| [«o]). 

Note that the state constraint boundary condition, i.e., the requirement that the trajectories stay 
inside the domain, is hidden in the control formula. We do not write it, however, explicitly, because, 
to the best of our knowledge, there is no general theory, as in [18], for state constraint problem with 
time varying and nonsmooth domains. Note that in our context we have no regularity properties for 
these domains. 

Given Cf +1 [uo] as in (J3T|) . it turns out that u* +1 [tio] is the minimal viscosity solution to 

uf +ljt [ U o] = \Duf +1 [u }\ 2 + R in Cf +1 [u ], 

(33) 
uf +1 [u ] = u in Cf +1 [« ] n (R d x {0}). 

Indeed using standard arguments from optimal control theory (see, for example, [H [2]), we may 
easily see that uf +1 [uq] satisfies the dynamic programming principle. The latter, as usual, implies that 
u i+i[^o] is a viscosity solution of (|33p . The proof of the fact that uf, i[i*o] is a minimal solution to 
(|3ll in Appendix [B 

The family {uf{uo]) i£ %+ $ >Q is nonincreasing in both i and S. Therefore there exists U [uo] > — oo, 
which is itself nonincreasing in 5, such that, as i — > +oo, uf [uq] \ U [uq] in M. d x [0,oo). 

Let U[uq] be the limit, as 5 — > 0, of the family (U [uq])s>q and, for \x > 0, consider the nonincreasing 
(in 5) family of sets 

n 5 [ Uo ] = p| nf[u ] and n[u - fj] = pi n 8 [u - ^}. (34) 

iez+ 5>o 

We have: 
Theorem 3.1. Let n £ be the solution to (CP, u e = eln(n e ) and assume ^). Then, for any \x > 0, 

u < U[uq] in R x [0, oo) and U[uo — /j] + \x < u in VI\uq — fi\. (35) 

Before we present the proof we remark that, by definition, u^[uo] = — oo in (C 4 - [«o]) • Therefore 
U 5 [u ] = -oo in (n 5 [u ]) c and, finally, U[u ] = -oo in (n[u ]) c = (f\ 5 CfM) C = (n*«%]) C , and, 
hence, 

u = — oo in (f][no]) c . 
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Moreover, since uf[-] > u m — 5 in Q, [•], by passing to the limit i — > oo and 5 — > we also obtain 

[/[•] > u m in «[•]. 

An important question is whether, as ^ — ?■ 0, C/[uo — //]—>• C/[tto]- This is, in general, not true. A 
counterexample can be found for u° = u m and R > 0. Then fq[tio — /^] cannot touch M. d x {0} and 
uf [no — u] = — oo. Therefore U[uq — u] = — oo for any /i > 0. On the other hand, u\ [uo] > u m and 
U[uo}=u\ 

We continue with the 



Proof of Theorem \3.1[ First we show by induction that, for all 6 > and i £ Z + , u < -u^fuo]. 

Since n* is a supersolution to (pQ), it follows from the comparison principle that n £ < n\ and, hence, 
u < u\[uo] = u 1 . 

Next we assume that u < ufluo], and, arguing by contradiction, we show, following an argument 
similar to that in Section 2, that u < uf +1 [uo] = — oo in (Q,f[uo]) c . 

To this end, suppose that, for some (xo,£o) G (^f[ n o]) C ) ^(^Oj^o) > — °°- Since u is upper semi- 
continuous, there exists a family (4> a )a>o of smooth functions such that u — <p a attains a strict local 
maximum at (x a ,t a ) and, as a — > 0, (x a ,t a ) — > (xo,to) 5 u(x a ,t a ) > u(xo,to), and, consequently, 
u(x a ,t a ) — > u(xo,to). It follows that there exist points (x ai£ ,t Qj£ ) £ (M rf x (0, oo)) where u £ — 4> a 
attains a local maximum and, as e — > 0, (x ai£ ,i Qi£ ) — > (xa.,£ Q ). 

Moreover, in view of (JH), at (x Qj£ ,t aj£ ), 

4> att -eA(f> a - \D<p a \ 2 - R< -exp((2e) _1 (u m -'u e )). 
Letting e — > yields, at (x Q ,£ a ), 

Q)t - |.D^ Q | 2 - -R < limsup(-exp[(2e)" 1 (n m - u £ (x a7£ ,t a>£ ))]). 

Since, by the definition of u, we have limsup e ^ u £ (x at£ ,t aj£ ) < u(x a ,t a ), the induction hypothesis 
yields that, for a small enough, u(x a ,t a ) < uf[uo](x a ,t a ) < u m — 5/2. 
It follows that 

limsup(-exp[(2e)~ 1 (n m - u £ (x a>£ ,t at£ ))]) = -oo, 

and, hence, at (x a ,t a ), 

4> a ,t - \D<j) a \ 2 - R < -oo, 

which, of course, is not possible because <f> a is a smooth function. 
Hence we have u = — oo in (£lf) c and, in particular, u = — oo on cftqfuo]- 

Next we show that 

u < u 5 i+1 [u ] = -oo in (Cf +1 [u ]) c . 

To this end, let (x,t) £ (Cf. , i[^o]) c \ (O^[uo]) c . Note that the existence of such a point means that 
(x,t) cannot be connected to M x {0} by a C 1 -trajectory staying in il^[uo]. Hence (x,t) belongs to 
a connected component V of wf[«o] = {{y,s) £ ^f[^o] '■ s <t}, which does not touch M. d x {0}. (See 
Figured!) 

Therefore d p V C <9f^[uo], where d p T> = {(y, s) £ dT> : s < i} is the parabolic boundary of V. From 
the previous argument we obtain 

u = — oo on d p V. (36) 
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As in (f2Tj) . for A > 0, we define w A by n e + exp (— ) = exp (^-)- Arguing as in the previous 
section, we deduce that, for all A > 0, 

uJ A = max(-A, u) and mintur 4 + A, w A - \Dw A \ 2 - R) < 0, 

and, in view of ([36]) . 

(min(uT A + A, w A - \Dw A \ 2 - R) < in £>, 
juT 4 = -A in <9 P £>, 

which admits —A + Cii as a supersolution for some Ci > 0. 
It follows from the comparison principle that, for all A > 0, 

u < w A < -A + dt in P. 

Letting A — > oo yields u = — oo in P and, consequently, u(x,t) = — oo. Observe that u = — oo in 
(Cf +1 [u ]) c implies that u = -oo on <9Cf +1 [u ] D (M d x [0,oo)). 

Finally we show that 

u<u 5 i+1 [u ] in Cf +1 [u ]. 

To this end, define z £ by n £ + exp ( x+1 £ - ) = exp (^) and notice that 

j r l . i / ( ^-"wN ^ , A , , / f uj +l [uo} - u £ \ \ 

z £ = u° i+1 [u \ + e In exp +11 = u £ + e In exp — ^ I + 1 I . 

It follows that 



z = max(«, Uj +1 [«o])- 



We claim that z is a subsolution of 



It - \Dz\ z - R < in C? +1 [uo]. (37) 



Indeed itf +1 [ito] is a subsolution to ([57]) in C,f +1 [uo] by definition. Moreover if, for some (x,£) E 
Cfj_ 1 [«o]) u(:e, i) 7^ — oo, using ([3]) and the stability of viscosity subsolutions, we find that u satisfies 
the viscosity subsolution criteria for (|37p at (x, t) . Finally, since the maximum of two subsolutions is 
always a subsolution, we obtain that z is a subsolution of ([57]) . 

We proceed by noticing that, since, in view of the above, 

u = -oo on <9Cf +1 n (R d x (0, +oo)), 

it follows that 

z = uf +1 [n ] on dC- +1 [u ] 

and, hence, 

z<uf +1 [u ] on <9C,f +1 [u ]. 

Therefore, using again the comparison principle for ([33]) . we obtain 

z<u s i+l [u ] in C.f +1 [u ] 
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and we conclude that u < n^ +1 [no]. 

Finally, since for all 5 > and i € Z + , we have u < n^[n ] it follows that, for all <5 > 0, u < 
linij^oo uf [uq] = II s [uq] . After letting 5 — > we obtain 

u < lim U S [u ] = f/[n ] in R d x (0,oo), 

which concludes the proof of the first part of the claim. 

For the second part we need the following lemma which is essentially a result from [7] that we adapt 
to our context (see also [3]). Its proof is postponed to the end of this section. 

Lemma 3.2. For all i £ Z + the lower semicontinuous function vf = max(n^[uo — u] + 25, u), is a 
super solution of 

' 4 t - \Dvf | 2 - R > in nf[uo-fj], 



vf = u in {u - u > u m - 5} n (R d x {0}). 

Since uf +1 [uo — u] is a minimal solution of ([33]) in C$ +1 [uq — fi] C Qf[uo — p] with 
Ui + i[uo — p] = uq — p on M. d x {0} (see Appendix |B]) , it follows that 

«H-i [u - p] < vf - p in Cf +1 [uq - fi] , 

and, hence, 

uf + i[u - /j] + n < max^ufluo - fi] + 25, u) in Cf +1 [n -u]. 

Letting i — > oo yields 

[/ [uq — p] + u < max([/ [no — u] + 25, n) in £1 [no — p\. 

Choosing p > 25 we also get 

C/ 5 [n -u] + 2(5<[/ 5 [no-/i] + n in fi 5 [n -u], 

and, therefore, 

[/ [no — u] + n < n in SI [no — u] . 

Finally letting 5 — > we obtain 

[/[no — u] + n = lim U [no — u] + p < n in O[no — //]. 
<5-»o 

We conclude with the 



(38) 



D 



Proof of Lemma \3.SX The key idea of the proof comes from [7J and [3] and relies on the property that, 
for concave Hamiltonians, the maximum of two supersolutions is supersolution. Here we reprove this 
fact in the context of semicontinuous supersolutions in a space-time domain. 

To this end, fix i € Z + and (x,t) G Cf [uq — p]. Since C\ \uq — p] is an open set, there exists p > 
such that B p (x,t) S C 4 [no — p], where B p (x,t) denotes the open ball of radius p centered at (x,t). 
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For a > 0, we define 

i4' a {x,t)= inf {uf[u - fi](y,s) + (2a)~ 1 {\x-y\ 2 + \t-s\ 2 )} and uf a '^ = uf a * xp, 
(y,s)eBp(x,t) 

where xp is a standard smoothing mollifier. 

Since u i ,a is an inf-convolution of the continuous function uf (see [2] ) , it is locally Lipschitz contin- 
uous and semi-concave with semi-concavity constant 1/a. 

It follows that Mj' is a smooth semi-concave function with semi-concavity constant 1/a and 

liminf u-' a,/3 (y, s) = uf[u - (jt](y, a). 



Finally, using Jensen's inequality and the concavity of the Hamiltonian, we obtain that, for some 
> 0, it/ is a smooth and, hence, a classical supersolution to 

4f P - eAut a>f) - \Duf a 'P\ 2 -R* X p> ~Ka - e/a in B p (x, t). (39) 



To prove ()38f) we show that the smooth approximations v^ e of vf in B p (x,t) given by 

6,a,p _|_ nf 8,a,P,e 

n £ + exp (^ Z±l) = exp (^ ) (40) 



are almost supersolutions to (|38p for a, /? and e small. Notice that in (|40|) we use 25 instead of <5. 
Replacing n e by exp (-*-^ — ) — exp (-* — ) in ([I]) we get 



ito e - /3 e ^7 = (#' ^ " eAvf a ^ £ - \Dv} aAe \ 2 ) exp^" 1 ^^) 



it 



S,a,p A S,a,P i n <5,ct,/3i2\ r — 1/ S,a,0 



i.t 



eAut' ^-{Du^ ,M | z )exp(e L {u i ' '^ + 25)), 



and, in view of (00 



S,a,p.e a S,a,P,e i n <5,a,/?,ei2 



!/ 



eA<' a '" - |£>t£ a 'f -R* X p) exp(e- l (u°> a ' p + 25- v°> a ' p ' e )) 
+ {R* X p-R) exp(e- 1 (^ a ' /3 + 25- vf*' 15 )) + R- fl £ n l J 2 ex^-e^vf"^) 

Using that, in view of flU, exp(e- 1 (nf' Q,/3 + 25 - v\' a ' Pfi )) < 1, and ([39D we find 

8, a 



v it H —eAvi — \Dv i | — R > —Ka — e/a 



+ {R* X p-R) exp(e- l (uf a ^ + 25- vf' a ' )) - (3 £ n l J 2 exp^ 1 ^) 

Define 

v 5 i ,a,P (y,s) = liminf vf a,P ' £ (y,s). 

Letting e — >■ and using the stability of viscosity supersolutions we obtain 



hit - l-°«i I 



(41) 
lim inf \{R *xp~R) exp{e^{uf a ^ + 25- vf a ^)) - $ e n x J 2 expOrV'"^)] . 

£— >0_ 
(»>s)->(5,s) 
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Recalling that u\ + 25 > u m + 6 in 0|[«o — A*], we deduce that, as e, a, /3 — > 0, in B p (x, t), 

Pen]/ 2 exp(-e- 1 ^ e ) = ^^/ 2 (n £ + exp^ 1 ^))" 1 < (l/2)/3 e exp(-(2e)- 1 n^' /3 ) -> 0. (42) 

Moreover, as e, /3 — )• 0, we also have 
# * X B ~ R -> 0, exp^" 1 ^'"^ + 25 - vf Q,/3 )) < 1 and uf (jj, s) = liminf vf^iy, s). (43) 

a,/3— >0 

(i/,s)->(!/,s) 



Using (j4Tj) , (j42|) , (j43j) , and the stability of viscosity supersolutions we find 

vl t - \Dvf\ 2 -R>0 in B p {x, t). 

Since all the above hold for all (x, t) E ^[uq — /z], it follows that the lower semicontinuous function 
vf is a supersolution to 

vl t - \Dvf\ 2 - R > in Q,f [u - p], vf = u for {u - (i > u m - 5} n (R d x {0}). 

D 

4 Constant rate 

Here we assume that the rate is a constant, i.e., 

R(x) = R in R d , (44) 

and, in addition, setting O = {x G M d : tto(x) > u m }, we have 

O = {x £ M d : n (x) > u m }. (45) 

We have: 
Theorem 4.1. Assume ^44\ ) an d i4°V - Then 

\im u e (x,t) = U[uo](x,t) locally uniformly in [M. x [0, oo) ) \ {(x,£) I U[uo](x,t) = u m } , (46) 

I |2 

O[u ] = {(x,i)| sup{- |a: ~ y| +^t + u (|/)}>n m }, (47) 

yao 4t 

and 

( 'su P2/e o [- l -^-+Rt + uo(y)} if (x,t)€Q[uo], 

-oo otherwise. 



U[u ] 



We notice that, if R < 0, then one can obtain (|46p from (|17|) and the dynamic programming principle. 
We also remark that, in particular, Theorem 14.11 shows that the limit of the family (u £ ) £> q is not, in 
general, given by (|17|). We refer to Appendix [X] for an explicit example. 
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(t,x) 



(0,o) (0,c) (0,6) 

Figure 2: The case with i?(x) = R a positive constant, {x € R| uq(x) > u m } = (a,b) and mq(') > ^ 



on an interval [a,b\. Then $7 = U^i a M{(i,t) 



\x-d\< 
it 



+ Rt + Uo(d) > u m }, the optimal trajectories 



are straight lines and U(x, t) 
is attained. 



\x-cY 

u 



+ Rt + Uq(c), where c is a point where the maximum in ([4*8 



Proof of Theorem \4-l\ When the rate R is constant, after one iteration of (|30p . pip and (|32p we find 
both the set O [«o] and the function U [uq], since for all i > 1, j > 2 and 5 > 0, ^[wo] = Cf[uo] = 
Q s [u } and uf[u ] = U 5 [u \. 

Indeed, every optimal trajectory in C 2 is a straight line connecting a point in £l 2 t° a point in 
I = {x G M rf : Uo(cc) > u m — 5} and, hence, it is included in f^. This follows from the observation 
that 



<f)(x,t) 



4/ 



i?t + Uq(c) 



is concave in (x,t) and, therefore, all the optimal trajectories of the points in Q 2 are included in Q^. 
It follows that r^2 = C|, ^2 = ^3 and consequently f^ = ^3- By iteration we obtain, for all i > 2, 
O^ = Of = O 5 and u\ = uf = £/ 5 . 
Using (j32j) and (J301) we see that, for all i > 2, 



0%] = ft? [«,,]= {(M) = sup{- 



|x-y| 



and 



U 5 [u ] = uf[u ](x,t) 



yClS 

jsup yeI s {- 



4/ 



+ Rt + u (y)} > u m - 5}, 



it 



+ Rt + u 



I —00 otherwise. 
It is easy to verify that (f33|) holds, since, for all * > 2 and 5 > 0, 



(y)} if (x,t)en s [u ], 



'i.t 



u }-\Dut[u ]\ 2 -R = in Of[« ]=Cf +1 [u ]. 



Letting 5 ->■ in (gU) and ((50]) we obtain (|47])-(j48]). (See Figure (2) 
We also have 



fl[u - fj] = P| O 5 [u - /i] = {(x, £) : sup { 



<5>0 



yeJ^ 



[a? — a/i 

\t 



+ Rt + u (y)} >u m + n}, 



with 



J» = {(x,t) : Uq > u m + fi}. 



(49) 



(50) 
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It follows that 

I |2 

\x — v\ 

U At>0 £l[u - p] = {(x, t) : sup { — \-Rt + u (y)} > u m }, (51) 

y&o 4t 



and 

\y-y\ 



lim U\uq — u] = < 



su P y eo {- 4* + fii + uo(v)} for (M) 6 U M>0 ^[ito - M]> 
— oo otherwise. 



(52) 



We also notice that 



I |2 I 12 

sup!- ' 3 ""/ 1 + fit + u (y)} = sup {- |X ~ y| + fit + u (y)}. (53) 

yeo 4* y6 o 4t 

Comparing gT]), (gSJ) with ([5T]). (|52]| and using ([53]) we deduce that 

lim U[u - p] (x,t) = U[uo](x,t) for U[uo](x,t) / u m , 

and, consequently, 

lim u £ (x,t) = U[uq](x, t) locally uniformly in (R x [0, oo) ) \ {(x, t) I U[uq](x, t) = u m } . 

D 



5 Strictly positive rate 

In this section we study the limiting behavior of the family (u e ) e> o when 

R > a > in M. d , (54) 

and show that, in general, the limit is not given by (|17p . 
For this we need to assume that, for sufficiently small p > S > 0, there exists ps :IM > such that 

lim lim ps a = and, if uo(y) > u m — 5, then sup uq(z) > u m — 5 + p. (55) 

<*-*°*->° ' lv-*l<Pd,„ 

Notice that it is important that p$^ is chosen independently of y. If uo G C 1 , ([55]) implies u m is 
never a local maximum of u. 
We have 

Theorem 5.1. Assume (54\ ) and 155]) , Then 

lim u £ = U[uq] locally uniformly in U At> o Q[uq — p\- (56) 

Recall that, in view of Theorem 13 .1\ we already know that lim e _>o u £ = — oo in f2[uo] c . 
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Proof of Theorem\in\ For h > h = £ + \\] £ + ^f, (x,£) G M d x [0, oo), i > 1 and ^ 5 > 0, we 
have 

u| [«o] (ar, *) < «f [«o - lA (x, t + h). (57) 

We postpone the proof of this inequality to Appendix ICl and we continue with the ongoing one. 
Letting i — > +oo and 5, [i — > we find, for all h > and £ > 0, 

U[u ](-,-)< Km U[u -ij](;-+h). (58) 

/i->0+ 

Hence, for all (x,t) £ U At >o^[no — fj], 

U[uo](x,t) < lim £/Tito — fA(x,t + /i) < u(x,t + /i) < u(x,t + /i), 

L r [no](x, £) < liminf u(x, t + h) < lim sup n(x, £ + /i). 
The definitions of « and u also imply that 

liminf u(x,t + /i) = u(x,t) and lim sup u(x,t + £t) = u(x,t). 

Combining all the above we obtain 

U[uo]<u<u in U M >o fi[i*o — A 4 ]- 

This last inequality and ([35]) yield u = u = U[uo] in U M >o ft[uQ — fi], and, hence, 

lim u e = U[uq\ in U, t>0 ^[n - //]. 

□ 



6 Conclusions 

We showed that the local uniform limit, as e — > 0, for the parabolic problem ([T]) with finite time 
extinction is naturally analyzed using the Hopf-Cole change of variables (|3|). The formal limit is 
the variant ([7]) of the standard eikonal equation. The new feature is the resulting quasi-variational 
inequality with an obstacle that depends on the solution itself. 

The quasi-variational inequality admits many solutions (see Appendix |A|) and the difficulty is to 
select the correct additional information. This is easy when the rate R is negative, as shown in Section 
[2j Indeed, in this case it is enough to enforce the Dirichlet boundary condition on the boundary of 
the unknown open set f2 where the liminf of the family (u e ) e> o is finite. This is due to the fact that, 
for concave Hamiltonians, the supremum of two supersolutions is still a supersolution. 

When the rate R is positive we do not have easy supersolutions at hand, and the answer is more 
elaborate. It requires an induction argument which allows us to identify again the limit of the family 
(^e)e>o- The key ingredient is a free boundary problem defined through the level set of the solution. 
The boundary condition for the resulting equation involves state constraints which leads us to study 
the problem using the related control problem. 

If the growth/death rate R changes sign, we can only bound from above and below the half-relaxed 
limits of the family (u E ) £> o by sub and supersolutions u and u respectively. 
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In terms of the biological motivation, our results qualitatively mean that the specific form of the 
survival threshold (a square root here) is irrelevant for the asymptotic problem. It also shows that the 
solution is deeply influenced by the survival threshold except when R is nonpositive. This confirms 
earlier numerical simulations in [14]. 

We conjecture that these upper and lower solutions are in fact equal and the correct setting (implying 
uniqueness) is to find a pair (u, Q) for which we can impose both Dirichlet and state constraints 
boundary conditions. Both establishing directly these boundary conditions for the half-limits of the 
family (u e ) e> o as well as developing a theory of state constraints boundary conditions for time varying, 
non-smooth domains are challenging mathematical issues. 

A Non-uniqueness 

To explain the difficulty associated with ([7]), we present here counter-examples for uniqueness and 
elaborate further conditions. Recall that the problem is to find pairs (u, £1) such that u is a viscosity 
solutions to 0. 

A first source for non-uniqueness is the value of u on dQ. Indeed assume that R and uq are such 
that there exists a unique viscosity solution u l of (|13p or, more generally, with u 1 defined in ([12]) and 
pip . For all rj > u m , we introduce the pair (u^f^) given by 

{n 1 (x,i) if (x,t) G il v , 
— oo otherwise. 

It can be easily verified that (w n , f^) is a viscosity solution of (|7|). In order to avoid this artefact, one 
can add the Dirchlet boundary condition ([S]) which appeared throughout our constructions. However 
in the next example we see that this Dirichlet condition is not enough to obtain uniqueness. In fact 
a state constraint boundary condition is hidden behind the property u 1 = — oo in the complement of 
Q v and we do not take it into account here. 

Let 

R(x) = 1 and Uq(x) = —x . 

A simple computation shows that the solution u 1 to (|13|) is given by 

u (x, t) = t 
Therefore the first truncation of u , given by 



x- 



1 + 4* 



u(x,i) 



b l+At iul L 1+44 — u *< 

— oo otherwise, 



with 

= {(x,t) : u(x,t) > — oo}, 

is a viscosity solution of ([7[). As a matter of fact this is the maximal subsolution to (]7|), (|5|) but it does 
not satisfy the state constraint boundary condition. To see this choose u m = —0.04. The point (1,2) 
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is included in Q since n(l,2) = 0.2 > —0.04. The optimal trajectory associated to this point, giving 
the value 5(1,2) = 0.2, is the straight line connecting (0,0.4) to (1,2). But n (0.4) = -0.16 < -0.04. 
So the point (0, 0.4) is not included in £1. Therefore a part of the optimal trajectory of the point (1, 2) 
is not included in £1. Hence 5 does not satisfy the state constraint condition. 

Following the arguments in Section U] we can find a viscosity solution to ([7]) and ([8j) . Indeed using 
(|4"8|) it is possible to compute explicitly the function U[uq] = lim^o U [uo] = lim^o n^fno] to find 



u(x,t) 

k — oo otherwise, 

with 

Cl = {(x,t) : u(x,t) > — oo}. 

From Theorem 14.11 we know that 5 is indeed the pointwise limit of the family (u e ) e >o outside the 
exceptional set {(x, t) : u(x,t) = u m }. 

However, in general 5 ^ u. Consider, for instance, the value u m = —0.04. Then 

5(2, 1) = 0.2, 5(2, 1) =0.15, 5(2.21,1) =0.02, 5(2.21, 1) = -oo. 

and, consequently, Cl £ fi. 

On the other hand, according to Section HI the state constraint boundary condition is satisfied for 
5, which motivates our conjecture in Section [6j 

B ufluo] is a minimal solution of (1331) in Cf[«o] 

Here we prove that nf[tto] is a minimal solution of (|33j) in Cf[no] by considering a supersolution 
w G Cf [no] of ([33]) and showing that 

ufluo] < w in C,f [n ]. (59) 

To this end, we fix (x,t) G Cf [uq\ and assume that (7(-), •) : [0, t] — > O^_ 1 [no] is a C 1 -trajectory with 
(j(t),t) = (x, t). Since Cf [uq\ is the set of points that can be connected by a C 1 -trajectory in £2f_ 1 ['Uo] 
to some point in M. d x {0}, it follows that 7 is included in Cf [uq\. 

For the supersolution w, we define, for s G [0,t], the (clearly) lower semicontinuous function f(s) = 
w('y(s),s) and we observe that (p is a viscosity supersolution of 

y/>-bl +jR ( 7 ) in ( ,t). (60) 

We postpone the proof of this claim to the end of the present paragraph and we proceed noticing 
that the function 

m = f ( - ^^ + R{l{s)))ds + «,( 7 (0), 0), 
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is a subsolution of ()60p . Then using the standard comparison principle of viscosity solutions we obtain 
w(x, t) = <p(t) >J\~ lj ^~ + R(l(s)))ds + n ( 7 (0)), 

and, since this is true for any C 1 -trajectory 7 and any (x,t) £ C\ [uq], ([59]) follows. 

It remains to prove (|60p . Let £ C 1 ((0,t)) be a test function, assume that i is a strict minimum 
point of ip — (j> and consider the function 

Fpiy, t) = w(y, t) - cf>(t) + |y ~ 7 2 (t)|2 + (t - t) 2 , 

A* 

which attains a local minimum at a point (y^,^) such that, as [i — > 0, 

i t^O and l ^~ 7 fr )|2 -X). (61) 

Since w is a supersolution to (|33p . we have 

^u4Mlli viuort ^ ^ | 2(^-7(t M )) |2 

<M*/JH o 7w + 2 (i M -i) ^ — — 2 — + R {y»)- 

It is immediate that 

0%) + 2(t„ - t) > -LTML 2 + fl(l , ), 



and, after letting // — >• 0, we conclude using ([6T 

C The proof of (1571) 

We prove by induction on i that, for all /i > h = ^ + \\J ' £ + ^, i > 1, S > 0, and (x, t) £R d x [0, 00), 

uf[uo](x,t) < uf[uo - fj](x,t + h). 

Recall that itf [no] = ^ 1 [«o] and u\[uq — /j] = tt^tto — /1] = u^uo] — /J,, where u^uo] is the solution of 
(fT3l) . Moreover (J53J) yields 



u [uq]{- , t) + ah — [i < u [uo](-,t + h) — \i = u [uq — fi](-,t + h). 
Therefore, for all h > h > /u/a, we have 



and, consequently, 
If, for all h > h and t > 0, 

it follows that, for all h > h, 



UiM(-,t) <nf[u -/i](-,t + /i). 

uf[«o](-,t) <ui[uo-n](-,t + h), 

nf[uo] + he t cnf[uo-fi], (62) 
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where e± is the unit vector in the direction of time axis. 

Fix (x,t) G Cj- +1 [uo] C ^f[^o] and let 7 be a C 1 -trajectory in f2?k\to] connecting (x,t) to a point 
(y, 0) with uo(y) > u m — 5. It follows from ([55]) that there exists z G M. d such that |z — y\ < p$^ and 
uo(z) > u m — S + /i. Without loss of generality we can take 'Uo(z) > uo(y)- 

The claim is that the trajectory 7 : [0, t + /i] — )■ R defined by 

{h~ 1 s{y-z) + z if 0<s</i, 
(63) 
7(5-/1) for h < s <t + h, 

is included in r^[«o — /x]. Indeed notice that the choice of h yields, for all h > h, 

\v - A 

+ a/i > n > 0. 



4/; 

Consequently, it follows from (|54p and the choice of z that the straight line connecting (y, /i) to 
(z, 0) is included in £1 [uq — fj] = r\jfij[uo — [i], and, in particular, in Qf[uo — //]. Therefore, for all 
< s < h, the point (j(s), s) is included in Qf [uq — fi]. 

Moreover using (|6"2j) we find that (7(5), s + /i) G f2f[iio — fi] for all s > 0. Hence, for all h < s, 
(7(5), s) G Sq [uq — ju], and we conclude that 7 is included in fi^[«o — A 4 ]- 



Next write 



( _ m£^ + m{s)))ds + Uq(z) _ M = / (_M. + ij( 7(s ))) ds (64) 

4 Jo 4 



+ j\-h^l + l?(7(5)))d a + u (z) - p. 



It follows that 

r \ l7(«)| 2 



y (-^^- + R(l{s)))ds + u (z) -fi> u (y). 
If this is true, then using (I30p . (|64p and (I65p we obtain, for all h > At and A > 0, 

^ + iH(',i) < u s i+1 [u - /j](-,t + h), 
and we deduce (1571). 



It remains to prove (J65J). Since i? > a, in view of (|63p . we have 

f h |7(s)| 2 _ |y-^| 2 
wo(^)>«o(y), / ( j \- R(l{s)))ds + u (z) - n > — \- ah + uq(z) - fi, 

and, for all h > h, 

\y — A r 

+ ah > fj,. 



Ah 
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